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Abstract: We discuss a diffusively perturbed predator-prey system. Freedman and 
Wolkowicz sliowed tliat tlie corresponding ODE can liave a periodic solution that 
bifurcates from a homoclinic loop. When the diffusion coefficients are large, this so- 
lution represents a stable, spatially homogeneous time-periodic solution of the PDE. 
We show that when the diffusion coefficients become small, the spatially homogeneous 
periodic solution becomes unstable and bifurcates into spatially nonhomogeneous pe- 
riodic solutions. The nature of the bifurcation is determined by the twistedness of an 
equilibrium/homoclinic bifurcation that occurs as the diffusion coefficients decrease. 
In the nontwisted case two spatially nonhomogeneous simple periodic solutions of 
equal period are generated, while in the twisted case a unique spatially nonhomoge- 
neous double periodic solution is generated through period-doubling. 

Key Words: Reaction-diffusion equations; predator-prey systems; homoclinic bi- 
furcations; periodic solutions. 
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1. Introduction 

Suppose that the ODE system 

U' = F{U,k), U eR"^, keR, F eC^{R^ xR;R^), (1.1) 

has a homochnic solution U = q{t) when the parameter k = k^. Assume also that 
for fcoo — e < ^ < ^oo, there is a stable periodic solution U = p(t, k) bifurcating from 
q{t). We study the diffusively perturbed system 

Ut = DU^^ + F{U,k), 0<e<l, 

(1.2) 

[/5(t,0) = [/5(t,l) = 0, 

where D = diag{(ii, (^2} is a positive diagonal matrix. The boundary conditions 
ensure that U(t,^) = q(t) or U(t,^) = p(t,k) is still a solution for system (|1.2D. 
Results from indicate that when the diffusion coefficients are large, these 

spatially homogeneous solutions (to be called SH solutions in the sequel) are stable. 
However, when the diffusion coefficients become small, SH solutions may lose stability 
and bifurcate into spatially nonhomogeneous solutions (to be called SN solutions). 

Such a bifurcation can create spatially nonhomogeneous patterns. Existing litera- 
ture on pattern generation concentrates on small patterns generated through bifurca- 
tions of equilibria, or traveling waves constructed using transition layers p5|, 0. The 



mechanism of pattern generation studied in this paper is fundamentally different. 

Since many ODE models are approximations to more realistic models where diffu- 
sion is present, examples that lead to systems (|1.1| ) and ( |1.2|) are plentiful. Freedman 



and Wolkowicz, ||T^, ^ , studied a two species predator-prey system that models the 
group defense of prey against predation. They found a homoclinic solution q{t) at 
a certain parameter value k = k^. The homoclinic solution bifurcates into a long 
period solution p(t, k) when k^o — e < k < koo- Suppose now diffusion is added to 
the system as in ( |1.2| ). The region of (di, 0^2) where pit, k) is stable has been studied 
in [|18], but the bifurcation when parameters cross the boundary F of the region has 
not been discussed. The purpose of the present paper is to discuss the bifurcation of 
p(t,k) when (^1,^2) crosses F. The bifurcation of q(t) into SN homoclinic solutions 
will be presented as the limit when the period is infinity. 

The creation of SN periodic solutions is caused jointly by the homoclinic bifurcation 
in (|1 . Ij ) and an equilibrium bifurcation in PDE modes in (|1.2|) . Let k = k^o such that 



Ll| ) possesses a homoclinic solution q(t) asymptotic to a hyperbolic equilibrium E. 
It is easy to find a curve F in the {di, (i2)-plane where the linearization of (|1.2|) at E 
has a simple zero eigenvalue, and no other eigenvalue on the imaginary axis. When 
{di,d2) crosses F transversely, the equilibrium E of ( p..2|) loses the hyperbolicity and 
two SN equilibria bifurcate from it. To describe bifurcations of q{t) and p{t, k) when 
{di,d2) crosses F, we need the concept of the twistedness of the homoclinic solution 
q(t). Let (pc be a unit eigenvector corresponding to the zero eigenvalue, unique up to 
the multiplication by —1. It can be shown that the linearization of ( p..2| ) around q{t) 
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has a solution (j)(t) that approaches 0c as t ^ — oo, and approaches c*0c as t ^ +00. 
Here c* is a scalar function of (^1,^2)- The limit of the solution (f){t) as t —>■ —00 is 
in fact a tangent vector to Wi'^^^{E), transverse to the unstable eigenvector. See Fig. 
1.1. We say that the homoclinic solution of ( |1.2| ) is twisted if c* < 0, nontwisted if 
c* > and degenerate if c* = 0. 

All the three cases have been found in Freedman and Wolkowicz's example by the 
numerical computation of c*. See Figure 7.4. An equivalent definition of twistedness 
will be given in §3. I have recently found that the change of twistedness is a generic 
phenomenon when diffusions are added to an ODE system that possesses a stable 
homoclinic orbit. But the proof will require a separate paper. 



Fig. 1.1 

The bifurcation of p{t, k) is determined by the twistedness of the homoclinic solution 
q(t) at the point where (^1,^2) crosses F. Roughly speaking, two SN simple periodic 
solutions of equal periods are generated in the nontwisted case, while in the twisted 
case a unique SN symmetric double periodic solution is generated through period- 
doubling. A periodic solution U(t,^) is said to be a simple periodic solution if its 
trajectory in a function space stays near the orbit of q(t) and hits a cross section E to 
the orbit of q(t) precisely once. It is said to be a symmetric double periodic solution 
if it hits S precisely two times and satisfies a symmetry condition U{t + T, ^) = 
U(t,l — C,). Here 2T is the period of U(t,^). Finally, when q(t) is degenerate, it is 
possible to cross F in a way that no SN simple or symmetric double periodic solutions 
are generated. We do not discuss the existence of other types of SN solutions in 
this paper due to technical complications. The homoclinic twist bifurcation at a 
hyperbohc equilibrium was discovered in and later studied in flSl and But 
the homoclinic twist bifurcation discussed in this paper is new even in the ODE 
context. 



In a separate paper, II^, we show how our method can be used to prove the stability 
of the SN periodic solutions. 

System ( p..2|) will be studied in the intermediate spaces Da{0) and Da{0 + 1), 
< 6* < 1. These function spaces allow solutions of ( |1.2| ) to have so-called maximal 
regularity, and are normally used to study fully nonlinear parabolic equations . Our 
system is not fully nonlinear, but to prove the smooth dependence of the solutions 
on di and (i2, which are the coefficients of the highest derivatives in the equations, we 
need to use the maximal regularity of the solutions. 

Solutions of (|1.2|) satisfy a reflection symmetry about the midpoint of the domain 



[0,1], due to the special boundary conditions imposed there. For a function U{^) 
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defined on ^ G [0, 1], let (i?f/)(0 = ^(1-0, < ^ < 1. It can be verified if f/(t,0 is 
a solution to (|1.2|), so is RU{t,^). Consequently, if Ui is a SN periodic solution and is 
mutually disjoint with RUi, then we have a pair of SN periodic solutions related by 
the symmetry. On the other hand, if Ui has a nonempty intersection with RUi, then 
Ui is a 2T period SN solution satisfying U{t + T,C,) = RU{t,C,). The R symmetry is 
very important in this paper since we can show that local center manifolds and fiows 
on them respect the symmetry group R. We can also show that bifurcation functions 
derived by Lyapunov-Schmidt procedures are invariant with respect to R. A mapping 
/ : C[0, 1] ^[0, 1] is invariant with respect to R if f{RU) = Rf{U). 

Suppose now the Neumann boundary conditions in ( |1.2| ) is replaced by periodic 
boundary conditions. In addition to the refiection symmetry, there is also a rotation 
symmetry, i.e., U{t,^ + 6) is a solution if U(t,^) is a solution. The bifurcation picture 
is quite different. Spatially nonhomogeneous tori may be generated instead of periodic 
solutions. See [0. Periodic boundary conditions will not be pursued further in this 
paper. 

We introduce intermediate spaces Da{0) and Da{0 + 1) in §2. We then study 
invariant manifolds and their foliations in these spaces. These invariant manifolds 
and their foliations provide convenient coordinates to study dynamics of ( |1.2|) near a 
equilibrium E. Some important lemmas regarding the symmetry R are also presented 
there. 

The assumptions and the main results of this paper are given in §3. 

In §4 we prove some lemmas needed in the sequel. In §5, we use a Lyapunov- 
Schmidt type reduction to obtain a one-dimensional bifurcation equation whose solu- 
tions correspond to simple or symmetric double SN periodic solutions. Proofs of the 
main theorems are given in §6. In §7 we summarize our numerical results about the 



example from . 

Recently Sandstede |^ has constructed center manifolds around some homoclinic 
solutions. It is hoped that a center manifold that is tangent to q'{t) and (j){t) can 
be constructed some day. And it may bifurcate into a center manifold around the 
orbit of p{t,k). Thus the twistedness of the homoclinic orbit should be passed to 
the twistedness of the center manifold of the periodic orbit. And the bifurcation of 
the periodic solutions should be determined by an one-dimensional return map on 
this center manifold. Thus, we naturally expect to see the occurrence of simple or 
symmetric double periodic solution on this center manifold. However, I was unable 
to use the center manifold technique in this paper due to technical complications. On 
the contrary, the bifurcation function approach in this paper is easy to use. The trade- 
off is that only simple periodic solutions and symmetric double periodic solutions are 
discussed in this paper. Complete understanding of the dynamics near the homoclinic 
orbit is still open, especially around the degenerate point c* = 0. 
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2. Abstract parabolic equations, invariant manifolds and foliations 

The PDE system (|1.2|) is studied in the intermediate spaces Da{0) and Da{6 + 1). 
Let A be a densely defined sectorial operator that generates a Cq analytic semigroup 
e"^* in a Banach space X . For each < < 1, define Banach spaces 

DA^e) = {xeX\ W-mt^at^^'^Ae^'x = 0}, 
DA{e + l) = {xeDA\Axe DA{e)}. 

The norms are 

\\x\\0 = supo<t<i ||t^~^Ae^*x||A' + \\x\\x, 

\\x\\g+i = \\Ax\\g + ||x||;t'- 

Intermediate spaces DA{d + m), < 9 < 1, m ^ can be defined similarly. 

Throughout this paper, let X = [^^(0, 1)]^, A = (^^^^ ^ and Da = {u e 

[iJ^(0, l)]^,-ug(0) = ^^(l) = 0}. See for details about the intermediate spaces. 

Denote J-'{U,di,d2,k) = DU^^ + F{U,k). We can write (|1.2| ) as an abstract non- 
linear parabolic equation 

f/' = ^(f/,/i) (2.1) 



where /i = {di,d2,k) is a parameter. The solution for ( |2.1| ) with f/(0) = Uq will be 
denoted by U^{t, Uq). 

It is easy to see that : + m + 1) x ^ Da{9 + m), m > is C^. Also 

F : Da ^ X is C°° since F is C°°. The following existence theorem is from [pi. 



Theorem 2.1. For each Uq G Da{0 + 1), there exists r > so that \2. 1\ ) admits a 
unique solution U G C^{[Q,t]] Da{0)) H {[Q , t]; D a{0 + 1)). Moreover, U is a 
function of Uq and fi in the specified function spaces. 

Consider a time dependent linear system 

u' = A{t)u + f{t), 

u{s) = x, a<s<t<b, y ■ ) 



which comes from linearizing (2J.) around a particular solution. It is easy to verify 
that 

(1) for all t G [a, b], A(t) : DA{t) X is sectorial, DA{t) = F)a with equivalent norms; 

(2) for each < ^ < 1, 

DA^t){e + 1) = DA^r){e + 1) = DA{e + l), for all t, r G [a, 6], 
with equivalent norms; 

(3) A{-) G C([a, h];L{DA, X)) n C([a, h]-L{DA{d + 1), 
The following theorem is proved in [1]. 
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Theorem 2.2. Under the above conditions, there is a unique solution 

U e C{[s, b]; Da{9 + 1)) n C\[s, b];DAm 

to ( ^■^) for each x G Da{0 + 1) and f G C {[s , b]; D a{0)) . Denote the solution by 
U{t) = T[t,s)x when / = 0. Then T{t,s) extends to Da{0) by continuity. Finally, 
the variation of constants formula holds for solutions of ^2.^ ) with f ^ 0: 

U{t)=Tit,s)x + J'Tit,ofmc 

Using the variation of constants formula, many familiar results of ODE systems 



can be extended to ( |2.1| ), with almost identical proofs. The most useful ones in this 
paper are the smoothness of invariant manifolds and their foliations. 

After a shift of coordinates, assume that {0} G Da{0 + 1) is an equilibrium of ( |2.1| ). 
Let A = DuJ-^{0, fio) where D^ = Da and hq = {dio, d2o, koo)- Here (dio, d2o) G F so 
that zero is an eigenvalue for A. Let M = T — AU. System ( |2.1| ) can be written as 



U' = AU + M{U,ii). 



Let 



cr 



j[A) = a_ U o-Q U 0-+ 
Rea^ < —Am, 
Rea^ > \m, 
Re ctq = 0. 

for some Am > 0. Let X,Y and Z be the invariant subspaces corresponding to the 
spectral set a+, o"_ and ctq respectively. We will identify Da{6 + 1) with X x Y x Z 
by writing U = {x,y, z) ii U = x + y + z. Since X, Y, Z C Da{0 + 1), Rx, Ry and 
Rz are defined by restricting R to subsets of Da{6 + !)• It is also easy to verify that 
RX = X, RY = Y and RZ = Z. 

Since both Af : Da{0 + 1) x Da{0) and TV : Da x ^ A" are there 

exists a local center manifold that is for any u > 0, [pil, 13] . Using the method 



of 0, which treats semi-linear parabolic equations, we can prove that the center 
unstable and center stable manifolds are C^, and there exists invariant foliation 
of center unstable (center stable) manifolds by unstable (stable) fibers, if Lip Af is 
small. The smallness of Lip Af can be removed by modifying the equation outside a 
neighborhood of {0}, if we are only interested in local invariant manifolds and their 
foliations. In the following we show how to choose the modifier so that the reflection 
symmetry resulting from the Neumann boundary conditions will be preserved for the 
induced flow on the local center manifold. 

First we give coordinate-free definitions of global center unstable, center stable 
manifolds and their foliations. See for similar definitions given to semilinear sys- 
tems. 
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Definition 2.3. Let < Ai < A2 < Am- The global center stable manifold is defined 
by 

W^' = {Uo e Da^O + 1) I U^it, Uo) exists for alH > and 

\mt,Uo)\\e+i<Ce'^\ t>0}. 
The global center unstable manifold is defined by 

pyc" = {Uo e DA{e + 1) I f/,(t, Uq) exists for alH < and 

Define the global center manifold by 

For each Uq G W^' (or W""), the stable fiber W'{Uo) (or unstable fiber Ty"(t/o)) 
passing through Uq is defined by 

W'{U,) = {K)Giy^'^(O) I ||f/,(^,^/o)-f/*(^,Fo)||m<Ce-*^^ t>0}, 

w^{u,) = {K)GW^^"(o)| ||f/,(^,f/o)-t/,(^,Fo)||e+l<Ce*^^ t<o}. 

Obviously, W'^'^ is forward invariant and ly™ is backward invariant. is invariant. 
Also each point on (or ty^") belongs to one and only one stable (or unstable) 
fiber. The global foliations are forward or backward invariant in the sense that 

U,{t,W'{U^)) C W'{U,{t,U^)), t>0; 

f/,(t,iy"(f/o)) c iy"(f/,(t,f/o)), t<o. 

Let O C D^id + 1) be an open set containing the equilibrium {0}. Let T : 
Da{9 + 1) X M3 ^ Da{9) and : Da X M3 ^ A' be > 0. Assume that = 

in (9 X M^. Consider the system 

[/' = .F(t/,/i). (2.3) 
Definition 2.4. Assume that (|2.3|) has global invariant center stabel and center un- 



stable manifolds, and invariant foliations as defined in Definition p.3| . Local invariant 



manifolds VF^o", VF^o* ^^^d ^Zoc system ( p.l|) are the restrictions to O of the global 
invariant manifolds for system (|2.3| ). Local invariant foliations of Wi^^ and Wi^^ for 
system ( p.l|) are the restrictions to of the global invariant foliations of M^™ and 
W^'^" for system (O). 



Local invariant manifolds and local invariant foliations depend on the extension of 
JF to JF outside O and are thus not unique. Observe that Lip M is small inside O if 
the neighborhood O is small, due to the fact A/'(0) = and A/''(0) = 0. The purpose 
of extending JF to JF is to have a small Lipschitz number for M = T — AU outside O. 
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Observe that Da{0 + 1) C [if^(0, 1)]^ is a continuous injection and ||u||//2(o,i) : 
1/2(0, 1) \ {0} ^ M+ is C for any z/ > 0. Let ^ : K ^ M be such that 

ip{s) = 1 for |s| < 1 and = for |s| > 2, 
< ij{s) < 1. 

Let Af{U,fi) = Af{4'{\U\[H'2(o^i)]2/ p)U, fi), where p > 0. It can be verified that 

is C^, LipAf — s> as p ^ 0, and A/" = TV for ||f/||[//2(o_i)]2 < p. Recall that Da = {U G 
[H'^{0, 1)]^ : dxU = at X = 0, 1}. After checking the boundary conditions, we find 
that both M,N : Da^ Da are for any z/ > 0. Since Da{0 + 1) C C Da{9), 
Af : Da{0 + 1) ^ Da{0) is with Lip A/" ^ as p — > in such space. We can prove 
the following theorem by using the method employed in . 

Theorem 2.5. For any z/ > 0, there exists a small constant p > such that the global 
invariant manifolds for system ( \2.3i) are C embedded submanifolds in Da{0 + 1) if 
|p — Pol < P- Moreover 

W"' = {x = hi{y,z,p)}, 
ly- = {y = h2{x,z,p)}, 

where {x,y, z) & X x Y x Z and |p — po| < p. The function hi, i = 1,2, is in all 
the variables, with hi{0,0, po) = 0, Dhi{0,0, pq) = andDhi = 0{p). 
By a change of variable {x, y, z) {x^, y^, z^), 

= X -hi{y,z,p), 

y^ = y - h2{x,z,p), 

z^ = z, 

we can flatten these manifolds, 

= {x^ = 0}, 

= {yi = 0}, 

W" = {x^ =0,y^ = 0}. 

The change of variables preserves the symmetry, i.e., U = {x,y,z) {x^,y^,z^), 
implies RU = {Rx, Ry, Rz) {Rx^, Ry^, Rz^). 

Proof. The existence and smoothness of such hi, i = 1,2, can be proved similar to p. 
It can be verified that Af{Rx, Ry, Rz, p) = RJ\f{x, y, z, p). Thus RW = W, where W 
stands for l^"", or W. 

Let U = {x,y,z) e W''. Then RU e W"". Thus h2{Rx,Rz,p) = Ry = 
Rh2{x, z, p). Similarly, hi{Ry,Rz,p) = Rhi{y, z, p). It follows that 
RU {Rx\Ry\Rz^). □ 
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We will use the new coordinates {x^, y^, z^) to discuss invariant foliations for system 
(12:^ ). Let Uo e W. Then Uq = {0,0, zq) in the new coordinates. Denote W^Uq) 
and W^IUq) by W''{zq) and W^Izq). From Definition |2.3| , we can show that different 
points on do not belong to a same fiber W^Izq) or W'^{zo). Further more, we also 
know that all fibers on W^"^ (or W^^) have to intersect W^. 

Theorem 2.6. If p > is small enough, then the stable fibers W^{zo),zo G 
form an invariant foliation ofW'^'^ and the unstable fibers W'^{zo), zq G form an 
invariant foliation ofW^^, for \fi — fio\ < p- Moreover, 

W^izo) = {x^ = 0,z^ = zo + h3{y^,Zo,p) with h3{0,Zo,p) = 0}, 
W^{zo) = {y^ = 0, z^ = zq + hi(x^, zq, p) with h^i^O, zq, p) = 0}. 

The function hi, i = 3,4, is in all its variables, Dyh^iO, 0, po) = 0, D^h^iO, 0, po) = 
and Dhi = 0{p),i = 3,4. By a change of variables {x^,y^,z^) {x'^ , y'^ , z^) , 
which is defined implicitly by 

x^ = x^, 

y' = y\ 

z^ = z^ + h^{y'^,z^,p) + hi{x'^,z'^,p), 

we can flatten the fibers, so that 

W\z^) = {x^ = 0, z^ = zo], 
W^{z,) = {y' = 0, z' = zo}. 

The change of variables preserves the symmetry, i.e., if{x^,y^,z^) — > {x'^,y'^,z^) then 
{Rx\Ry\Rz^) {Rx^,Ry'^,Rz'^). 

The change of variable here does not affect the flow on W^. 

Proof. The existence and smoothness of hi, i = 3,4, are proved similar to that in ||^. 

i) Since W'{Rzo) = RW'{zo), and W^Rzo) = RWizo), we have h3{Ry\ Rzq, p 
= Rh^iy^ , Zq, p) and h^{Rx^ , Rzq, p) = Rh^{x^ , zq, p). It follows that {x^,y^,z^) 
(x^y^z^) implies that {Rx^ , Ry^ , Rz^) {Rx^ ,Ry'^ ,Rz^). 

ii) If = and y"^ = 0, then hz{0,z'^,p) = and hi{0,z'^,p) = 0. Therefore 
z^ = z^ on . The equation for the flow on is not changed. □ 

Define a function space A'n = {(m„ cos(n7r-), t>„ cos(n7r-)), (u„,t>„) G M^}. Obviously 
Xn is isomorphic to R^. Observe that n > 0} is dense in [-^^^(0, 1)]^. 

Recall that (|2.1|) comes from ( |1.2|) . The hypotheses on F will be specified in §3. In 
particular, they imply 

1) Z is one dimensional, spanned by an eigenvector in X\, corresponding to the 
eigenvalue A = 0; 

2) X is one dimensional, spanned by an eigenvector in Xq, corresponding to the 
eigenvalue A = A+. 
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We may identify Z and X with M. More precisely, let w be a unit vector in Z. For 
any 2; G there is a unique z G M such that z = zw. We will identify z with 1 and 
drop the over-bar. The same comment also applies to X. It can be verified that if 
X G X and z G Z, then Rx = x and Rz = —z. We use U ~ {x^, y'^, z^) to indicate U 
corresponds to [x"^ , y"^ , z^) in the new coordinates. 

Theorem 2.7. a) If U E Xq and if U ^ in the new coordinates, then 

= and G Aq. The converse is also true, 
b ) For system ( ^■^ ), the flow on has the form 

= 0, y^ = 0, 

where g{0, /i) = 0, D^igi^O, /io) = and gi^—z"^, jj) = —g{z^, jj) . 

Proof, a) In the original coordinates, U = x + y + z with x E X , y E Y and z E Z. 
HUE Xq, then it is obvious that z = 0, x E Xq and y G Xq. 

We first examine the changes of variables (x, y, z) {x^,y^, z^)) as in Theorem |2.5| . 
Consider the change of variable = y — h2{x,z,fi). When z = 0, the graph {z = 
0, y = /i2(x, 0,/x)} = W^^ n {z = 0} is one-dimensional. Now let us restrict the 
system to Xq, where E = {0} is hyperbolic. By the standard existence theorem of 
the unstable manifold for the ODE system, there exists a smooth function h such that 
W"^ = {y = h{x,fi)} for the restricted system. Clearly the graph {y = h{x,fi)} C 
{y = h2{x,0, fi)}. Since they are both one dimensional, we have that h{x,fi) = 
h2{x,0,fi). This proves that h2{x,0,^) G Xq. Recall that z = z^. Thus HUE Xq, 
z^ = and y^ E Xq, vice versa. 

We now consider the second change of variable {x^,y^,z^) {x'^,y'^, z"^) as in 
Theorem Since y"^ = y^, y"^ E Xq <^ y^ E Xq. If {0,y^,z^) is a point on 

W^{zo), then {0,Ry^,Rzi) is a point on W'^{Rzo). From Theorem |2.6| , compare 
the z^ coordinates, and observe that Rzq = —zq, we have R{zo + h^^y^ , zq, fi)) = 
-zq + hs{Ry^, -Zo,n). However, Rh^^ = -h^. Thus h^{Ry^, -^0,/^) = -hiv^ , zq, 1^)- 
Similarly, we can show h^^Rx^ , —Zq, fi) = —h^{x^,ZQ,ii). Therefore if y"^ = y^ E Xq, 
we have Ry^ = y^, h^{y'^,Q,ji) = 0, and hi{x'^,0,fi) = 0. In this case, we have 

= ^ z"^ = 0. 

By combining the two changes of variables, we have verified the assertions of a), 
b) The assertions = and y"^ = are obvious. If U(t) ~ (0, 0, z^{t)) is a solution 
on W^, so is RU{t) ~ (0, 0, -z\t)). Therefore g{-z^, fi) = ~g{z^, fi). □ 
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3. Assumptions and main results 



We assume that the ODE system satisfies the following hypotheses: 
Hi) F : M2 X M ^ M2 qoo^ 

H2) At k = koo, (|1.1| ) possesses a homoclinic solution U = q{t) asymptotic to an 

equilibrium E = E{koo). 
H3) At E{koo), the Jacobian matrix 

satisfies a + d > and ad — be < 0. 

Hypotheses H3) implies that E{k^) is hyperbolic with eigenvalues denoted by 
— A_ < < A+, satisfying A+ — A_ < 0. The equilibrium E = E{k) continues to 
exist for all k ^ k^o- We will suppress k if no confusion should arise. The homoclinic 
orbit is stable from inside since A+ — A_ < 0. Assume that the homoclinic orbit breaks 
in certain direction when k moves away from fcoo, so that periodic solutions bifurcate 
from q{t) for koc, — € < k < koo- More precisely, consider the linear variational equation 
of ( |Ll]) around U = q{t), 

U' = duF{q{t),koo)U, (3.1) 

and its adjoint equation 

^' = ~[duF{q{t),koo)]*^. (3.2) 
System (|3.2|) has a unique nontrivial bounded solution ^(t) up to multiplying by 



nonzero constants. It is known that \E'(t) ~ \E'oe~'^+* and q{—t) — E{koo) ~ (f)oe~^+^ as 
t +00 where \l'o (or (po) are left (or right) eigenvector of the matrix J corresponding 
to the eigenvalue A+. See [0]. For definiteness, assume 

lim ^(t)(g(-t) - E)e^'^+ = -1. (3.3) 

t — >H-oo 

We now assume that the breaking of the homoclinic solution q{t) is in the direction 
determined by 

H4) IZo ^ W ■ dkF{q{t), koo)dt > 0. 

From Silnikov (|3.3|) and H4) imply that there exists e > so that for k^o — e< 
k < koo, system ( |1 . 1| ) has a simple periodic solution p{t,k) which is orbitally near 
q{t) and is asymptotically stable. A more transparent relation indicating that the 
periodic solutions can only be found for k < koo with k — koo = 0(e~^'*'+) is given in 
[ P^ , where T is the period of p{t, k). There is a one-to-one correspondence between 
k and T. Moreover, there exists C > 1, independent of k, such that 

- dT - 



The proof of that can be obtained by the same method used in . 
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Consider eigenvalues for the linear variational equation around the equilibrium 
E(koo)- It can be verified that each eigenfunction must be in one Xn, n > with an 
eigenvalue A„ satisfying 

A + a + n^TT^rfi b ^ — 

c X + d + n^Tr^d2 J ~ 

The spaces Xn are defined in §2. For each Xn, denote the eigenvalues corresponding 
to the n-th Fourier mode by (A„i, A„2), with ReA„i > Re A„2- Based on a + d > 0, we 
have Re A„2 < 0. An n-th mode is unstable if and only if Re A„i > 0. The critical 
case A„i = occurs if 

(a + n'^n'^di){d + n'^n'^d2) = be. 

We can show that when decreasing {di, 1^2), the first mode loses stability before the 
other Fourier modes. (Theorem 3.1). Thus, we are interested in parameter values 
where An = 0. Define 

F = {{di, (^2) : (a + TT'^di){d + n'^d2) = be} 

Theorem 3.1. The first quadrant, M^, is divided by F into two regions: and Q- 
where {a + 7T^di){d + TT^d2) — be < and > respectively. 

(i) All > in Q^. If di and d2 are sufficiently small, then (c/i, 1^2) € 

(ii) Re All < in Q^. The region is unbounded. 

(iii) All =0 on F. 

(iv) Aoi = A+ > m M^. // (rfi, 4) G ^_ U F, then ReX^j < for (n, j) ^ (0, 1) or 
(1,1)- 

(v) VAii = (dd^Xu^dd^Xu) 7^ for (^1,^2) e F. In particular, dd^Xn < if 
d + Tx'^d2 > and dd^Xu < if a + Ti'^di > 0. 




Fig. 3.1 

Theorem p.l| will be proved in §6. Fig. 3.1 depicts F, and for all possible 
cases except for a possible permutation of di and d2. When be > 0, we may have 
a > 0, d > which is Case 2 in Figure 3.1. We may also have a > 0, d < which is 
Case 1 in Figure 3.1. It is impossible to have a < 0, d < since a + d > 0. The other 
possible case is a < 0, > which is obtained from Case 1 by symmetry. When 
be = 0, since ad — be < 0, we have ad < 0. Thus we have either a < 0, d > which is 
Case 3 in Figure 3.1, ora>0,(i<Oby symmetry. When be < 0, again ad — be < 
implies ad < 0. Case a < 0, > is in Case 4, the other case a > 0, < is obtained 
by symmetry. Observe in Case 4, when increasing d2, we can move from to Q^. It 
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is interesting to note that the equihbrium may become more unstable by increasing 
one of the diffusion coefficient. 

The following theorem was stated in p 



Theorem 3.2. For each positive (^1,(^2) € Q-, there exists a smooth function 
e*{di,d2) > such that the SH periodic solution p{t,k) is asymptotically stable in 
Da{0 + 1) ifkoo-e*<k<k^. 



Theorem |3.2| can be proved by using notions of exponential dichotomies and rough- 
ness of exponential dichotomies in Da{0 + !)• The proof is similar to the proof of 
Theorem 4.5 in Since those methods are quite different from the ones used in 
this paper, we will not give details here. 

The result in Theorem |3.2| is not very precise since e*{di, 6/2) — > as {di, ^2) F. 
For a given k (or period T), the loss of stability for p{t, k) does not happen exactly 
at F. 

To describe what happens near F, two new notions are introduced: 1) the stability 
of the equilibrium for the flow on Wf^^; 2) the twistedness of the homoclinic orbit 
when following q(t) from t = —00 to t = +00. 

When {di,d2) G F, An = 0, the equilibrium E{koo) has a one-dimensional center 
manifold Wi^^ that is tangent to the one-dimensional eigenspace corresponding to 



All = 0. The flow on Wi^^ is described in Theorem p77| . When An = 0, it has the 
form 

z' = —cz^ + h.o.t. 
X = 0,y = 0. 

H5) When All = 0, the equilibrium E is stable on Wf^^^E) in the sense that c > 0. 

Numerical computation in §7 shows that in Freedman and Wolkowicz's example 
the condition c > is valid for all {di, ^2) € F in the range specified by < vr^c/i < 3. 

Twistedness of the homoclinic solution q{t) has been described in §1. Because 
its importance, we will give a simple and equivalent definition. Let k = koo and 
(^1,^2) G F. Linearize ( |1.2| ) around q{t) we have 



U\t) = DU^^it) + duF{q{t), k^)U{t). (3.4) 



The subspace of the first Fourier mode Xi is invariant under ( |3.4]) . Since Xi is two 
dimensional, ( |3.4| ) on Xi reduces to an ODE on (ui,t>i). 

(3.5) 

where A{t) ^(00) = ( ^ ^ j f _^ ±00. According to a 

\ — c — vr U2 — cL I 

theorem in each solution of ( p.5| ) approaches a solution of the linear autonomous 
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equation, 

with an exponentially small error. Therefore, there is a unique solution {ui{t) , Vi{t)) 
to ( |3.5| ), up to multiplying by scalar constants, that approaches an eigenvector of the 
zero eigenvalue of ^(oo) as t — oo. Let U(t) = {ui(t) cos nx,vi(t) cos nx) be the 
unique solution of ( p.4|) that is in Xi and approaches an eigenvector 0c, corresponding 
to All = as t ^ — oo. By the same argument, when t +oo, U{t) approaches 
another eigenvector associated to An = 0, denoted by c*0c, where c* is a function of 
{di, G F. 

Definition 3.3. Let lim(_>_|_oo U{t) = c*(j)c- The homoclinic solution q{t) is said to be 
nontwisted if c* > 0, twisted if c* < 0, or degenerate if c* = 0. 



Remark. In Da{0 + 1), solutions of ( p.4| ) that approach 0^ as t — oo are not 
unique. They have the form U{t) = U{t) + Cq{t) where C is an arbitrary constant. 
Since q(t) exponentially as t — +oo, we have limj^oo U(t) = c*(l)c for any C G M. 
Therefore the twistedness defined in Definition is precisely the one given in §1. 



Let k = koo and {di, ^2) G F. From Theorem |3.1| , VAn 7^ 0. It is also obvious that 
VAii intersects F transversely at {di,d2). We can make a smooth change of variable 
B : {di, (^2) — > {i, m) in a neighborhood of F so that m = An and i is the arc length on 
F when m = 0, after assigning £ = to an arbitrary point on F. The new coordinates 
flatten F, i.e. F = {m = 0, £ G /}, where / C M is an open interval. 

For m ~ and i & I, we look for simple period T or symmetric double period 2T 
SN solutions, where T >t,t being a large constant. In the parameter space (T, i, m) 
we want to find regions where such SN solutions exist. 

For {di, (^2) G F, B{di, 1^2) = {^o, 0), and the twistedness c* = c*{£o) is a function of 

Throughout this paper, assume that the hypotheses Hi)-H^) are satisfied. 

Theorem 3.4. For each Iq G /, c*(£o) 7^ 0, there exist a large constant t > and 
an open set (D C containing {£0, 0), the size of which depends on io, such that two 
functions L{T, i, m) and r{T, i, m) can be defined for T >t and {i, m) G O, with 
the following properties: 

1) r(oo,4,0) = c*(£o), where r(oo,£o,0) = limf^+oo '"l^, ^0, 0); 

2) L{TJ,m) = e™^ + 0(e-"^), < m < a; 

3) £^{L{T,i,m)r{T,i,m)} > (or < 0) when c*{io) > (or<0). 

Moreover, the existence and uniqueness of simple or symmetric double SN periodic 
solutions to (2.1) are determined by the following conditions: 

i) If < c*{£o) < 1 or 1 < c*{io), then there is no simple period T SN solution 
when < L{T,£,m)r{T,i,m) < 1; there are precisely two simple period T SN 
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solutions Ui{t,^) and U2{t,^) when L{T,i,m)r{T,i,m) > 1. The two solutions 
are related by U2{t,^) = Ui{t, 1 — 0- 

ii) //c*(£o) = 1; then there exist two simple period T SN solutions when 
L{T,i,m)r{T,i,m) > 1. There exists 6 > such that the number of solutions 
are precisely two when L{T, i, m)r{T, i,m) > 1 + 6 for some 6 > 0; and there is 
no simple period T SN solution when < L{T, i, m)r{T, i,m) < 1 — 6. 

iii) If —1 < c*{io) < or c*(£o) < —1; then there is precisely one SN symmetric 
double period 2T solution U{t,C,) when L{T, £,m)r{T,i,m) < —1. There is no 
such SN period 2T solution when —1 < L{T,£,m)r{T,£,m) < 0. 

iv) If c*{£q) = — 1, then there is at least one symmetric double period 2T SN solution 
when L{T, i, m)r{T, i,m) < —1. Such solution is unique when L{T, i, m)r{T, i, m) 
< —1 — 6 for some 6 > 0. There is no such SN period 2T solution when 
-1 + 6 < L{T, i, m)r{T, i,m) <0. 

Corollary. When m > 0, there is a pair of SN equilibria Ei, E2 bifurcating from E. 
The results above also show the bifurcation of SN homoclinic or heteroclinic solutions 
asymptotic to Ei and/or E2 as a special cases when T = 00. If c* ^ 0, the curve 
Lr = 1 is identical to F. When crossing F at a point where c*(£o) > 0, the bifurcation 
of a pair of homoclinic solutions, each asmptotic to Ei or E2 occurs. When crossing 
F at a point where c*(£o) < 0, the bifurcation of a pair of heteroclinic solutions 
connecting Ei and E2 occurs. 

Theorem 3.5. For each io E I with c*(£o) = and ^c*(£o) 7^ 0, there exist constants 
e > and t > such that functions i*{m), \m\ < e and 6{T) = ce""*"^, T > t for 
some c > can be defined. If \i — £*{m)\ < 6{T), \m\ < e and T > t, then there 
is no simple period T or symmetric double period 2T SN solution to (2.1), inside a 
((5(T))^/^ neighborhood of the orbit of q{t). 



Theorem |3.4| provides fairly accurate information about the bifurcation to simple or 
symmetric double periodic SN solutions when crossing the curve -L(T, £, m)r(T, m) = 
1 not near the points c*{£q) = ±1 or c*(£o) = 0. First, Theorem 3) assures that 
L{T,i,m)r{T,i,m) is monotonic in term of m. From the asymptotic forms 1) and 
2), it is also clear that the sign of Lr — 1 changes when m is increased form nega- 
tive to positive, provided that T is large. . When crossing the curve Lr = 1 near 
c*(^o) = il) bifurcation to simple or symmetric double periodic SN solution will oc- 
cur but the precise moment is unknown. Our method does not predict the existence 



or uniqueness of such solutions in a narrow strip around Lr = 1. Theorem |3.5| , on 
the other hand, assures that when c*{io) = 0, we can pass m = through a small 
tubular neighborhood of i = i*{m) without creating any simple period T or symmet- 
ric double period 2T SN solution. The size of the tubular neighborhood shrinks to 
zero as T ^ +00. 

, . r^, it:— n , are 



The regions in the (^1,^2) plane mentioned in Theorem |3]^ and Theorem p3 
depicted in Figure 3.2 where we assume that L{T, i, m) = e"^'^, r{T, i, m) = c*{i) and 
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£*{m) = 0. In the shaded area, the existence and uniqueness of a simple period T (or 
symmetric double period 2T) SN solution is guaranteed except near c*{i) = ±1. The 
tubular neighborhood near £ = where crossing m = without causing bifurcation 
to simple or symmetric double period SN solution is also shown. A sketch of all kinds 
of homoclinic, heteroclinic and periodic solutions is in Figure 3.3. 



Fig. 3.2 



Fig. 3.3 
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4. Some lemmas 



The results in Lemma ^?T| is our major tool to study a solution U{t), < t < to, that 
stays in a small neighborhood of a nonhyperbolic equilibrium. Following an idea of 
Silnikov, we show that if to can be arbitrarily large, U{t) = {x{t),y{t), z{t)), < t < to 
is determined by, and depends continuously on its boundary values: y{0),z{0) and 
a;(to). Using exponential dichotomies we can easily show x(t) = 0{e~°'^^°~^^) and 
y{t) = 0{e~°'*) for some a > 0. However, in the center direction, the flow is not 
exponentially decaying either moving forward or backward. Following the approach 
of [0, we will compare the Z coordinates of U{t) with a (nonunique) solution f/o(t) 
on Wig^{E). Let P^, Py, be the spectral projections from Da{0 + 1) onto X, Y, Z. 
In the fiat coordinates, we show that \\Pz{U{t) — Uo{t))\\ is small and approaches zero 
uniformly for t G [0,to] as to +oo. If we are interested only in dynamics in the Z 
direction, U{t) can be replaced by Uoit) on Wi^^{E) with a very small error. 

It is also clear that the smallness of PziUoit) — U{t)) strongly depends on a good 
choice of coordinates. Since a; (to) and y{Q) are not small as to — > cxd, a undesired 
change of variables may destroy the smallness of Pj(f/o(t) — U{t)). 

Let O C 'Dji^lO + oo) be a small neighborhood of an equilibrium U = where 
the fiat coordinates introduced in 52 are used in O. We now consider the abstract 



parabolic equation (|2.1|) written in the flat coordinates, 

x' = Aix + gi{x,y,z,n), 

y' = A2y + g2{x,y,z,fj,), (4.1) 
z' = A3Z + gs{x,y,z,ij,). 

Here Ai = A\x, A2 = A\y and A3 = A\z. Rea{Ai) > Am > 0, Rea{A2) < -\m < 
and KealA^) = 0. The functions gi,i = 1,2,3, are C^, > 2 in all the variables. 
Since the coordinates are fiat, it can be verified that 5^1(0, y, z, fi) = 0, g2{x, 0, z, fi) = 
and 5(3(0, y, z, /i) = g^ix, 0, z, fx) = 53(0, 0, z, /i). Moreover, Dugi{0, 0, 0, ixq) = 0, i = 
1, 2, 3. The equation for the flow on the center manifold is 

z' = A3Z + g3iO,0,z,fx). (4.2) 

Let <l>(t,zo,/i) be the solution map for (|4.2| ), with $(0,2;o,/i) = zq. We have the 
following 

Lemma 4.1. For any uq, P > with < P < ao < \m, there exist positive constants 
^M, ^M, f^M tm with the following properties. The constant is small enough so 
that {U = {x,y, z) \ \\x\\x < ejv/, WvWy < cm, \\z\\z < cm} C O. If < /iM,to > t^ 
and zq G Z , satisfying 



\\^t,ZQ,fi)\\z<eM /ortG[0,to], 
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and if \xq\ + \yo\ < 6m , then there exists a unique solution U{t) E O, t E [0,to], to 
equation ( \4.1\ ), satisfying the boundary conditions 

x(to) = xq, y{0) = yo and z{0) = zq. 

The solution can he written in the form 

U{t) = (x^t), y%t), + z^t)), < t < to, 

where = ^{t,Zo,n), z^(0) = 0. w^{t) = w^{t;to, Xo,yo, zq, fx), w = x,y or z, 
are functions in all the variables if Qi, i = 1,2,3 is . Moreover, let r he a 
multi-index with < |r| <u — 1. Suppose ai satisfies Q < (3 < ai < — \r\(3 . Then 

\D''x^{t)\x < Ce"i(*-*o), 
\D'-y\t)\Y < Ce-'^i*, 



The proof for Lemma [4.1| in the ODE case can be found in |T^, |T^. The proof 
for systems of abstract parabohc equations is similar and will not be rendered here. 

Since the small eigenvalue An = m and since the flow on the center manifold is 
odd, we can rewrite ([4.2|) as the following, 



z 



mz — cz^ + z^hi{z, /i). 



(4.3) 



Here c = c(yu) > due to H^), and < hm- The function hi is for all z/ > and 
hi{—z,fi) = hi{z,n). Equation ( [4.3|) has three equilibria z = and z = ±ze, where 

Ze ~ Jm/c provided that m > and m is small. 



In Lemma [4.2| and Lemma |4.3| we present some estimates on the function 
(^{t, Zq, ^) / Zq which measures the degree of expansion or contraction on the center 
manifold. The importance of these estimates will be clear in the next two sections 
where bifurcation functions and their approximations are introduced. The proofs are 
technical and can be skipped on the first reading. In fact, the results in Lemma O 



and Lemma are easy to verify for the truncated equation 



mz — cz 



All we try to show in these lemmas is that the perturbation term z^hi{z, fi) does not 
change the solution significantly. 

Let e > be a small constant. By plotting the phase diagram of ([4.3|) on (— e, e), 
see Figure 4.1, it can be verified that \^{t, zq, fi)\ < e provided \zo\ < e, and m and 



/iM are small. In Lemma [4^, we show ^{t, zq, ^)/ zq is monotonic with respect to zq 
in (0, e) if t > is fixed. We also give formulas that will provide some lower bounds 



on 



in the future. 



Fig. 4.1 
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Lemma 4.2. There exists e > such that 

if < \zo\ < e. Moreover, we can show the following. 
(i) lfm<0, then 



-signzQ 



_d_ 

dzo 



Ciizi-'^' 



m - czl + z^hi{zQ, n) zl' 

where Ci is a function of zq. Ci ^ c if m and e are small. 
(ii) If m > and Zq ^ z%, then we have 



_d_ 

dzo 



C2{z'o-<l>' 
Ze ~ Zq 



2 ' 
Zo 



where C2 is a function of Zq. (72 ~ 1 if m and e are small. 
(iii) If m > and z^ = z"^, then 



_d_ 

dzo 



Zo 



-m't 



where —m' = 4^ [mz — cz^ + z^hi{z, /i)] 



z=ze 



l)/^0, 

—2m if m and e are small. 



Proof. Since $(t, zq,/^) is an odd function of zo, it suffices to consider zq > 0. Let 
z = ^{t,zo,n). It is well known that U{t) 



■i^ satisfies the linear variational 



equation for ( |4.3| ) with U{0) = 1, so does . Therefore, they must be 



identical. Using ( [4.3| ) to replace "'^''g^"'^-' , we have 



d^{t,zo,fj.) 

at ■ 

9$ mz — cz^ + z^hi (z, jj) 
dzQ mzo - czq + z^hi {zq, /i) ' 



d /$ 



dzo \zo 



mz — cz^ + z^hi {z, /i) Zo 



- 1 



Zo 



[ mzQ - czq + ZqHi^Zo, jj) z 
c{zl - z'^) + z'^hi{z, /i) - Zo^i(zo, /i) 



(4.4) 
(4.5) 



^2 
Zq 



m - czq + ZqHiIzq, fi) 

Since hi{z,fi) is an even function of z, we have 

z'^hi{z,n) - z^hi{zo,fi) = CsizQ - z^), 

where C3 is a function of zq and is small if both z and zq are small. This proves (i). 
For any fixed t > 0, m > 0, let zq ze, then 2; = $ ^ ze- From (4.4), and the 

Ze, t fixed, we have 

mz - cz^ + z^hi(z, /i) 

hm o ; , = e , 

^o^2B mzo - czq + 2;o/ii(zo, /«) 



fact — 



e as 2:0 
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where —m' = §^[mz — cz^ + z^hi{z, fi)]\z=ZE ~ —2m, since ze ~ \Jm/c. Therefore, 
(iii) follows from the first line of (4.5). 

When m > 0, since ze is a nonzero equilibrium, m — cz'^ + z^hi{zE, ft) = 0. 
Therefore 

m- cZq + ZQhi{zo,^) 
= m- czl + z^hi{zo, /i) - [m - cz% + z%hi{zE, /u)] 
= Cii^Z'^ — ZQ), 

where C4 is a function of Zq and is close to c. From this, (ii) follows from the second 
line of (4.5). 

When m < 0, zero is an attractor on Wf^^. If zq > 0, then Zq > Thus 
^[J] < 0, based on (i). When m > 0, zq = ze, ] < 0, based on (iii). When 

m > 0, Zq Ze, Zq > 0, Ze attracts zq. From the phase diagram, see Figure 4.1, it is 
clear that Zq — $^ and z'^ — Zq always have different signs. Thus ^[^] < 0, based 
on (ii). □ 

In the next lemma, we derive some estimates on the rate of contraction or repelling 
for the equilibria zero and/or ±ze on the center manifold. These estimates are to be 



used in conjunction with Lemma 4.2 



Lemma 4.3. There exist 'z > and fn > with the following properties. 

(a) Let \zq\ < z, < m < fn and $ = $(^0, ZQ,fi). For any tq > 0, if mto > tq > 

then there exists rj = tjItq) > such that either 

zl < (1 - V)^^ or |$2 - 41 < (1 - r^)\zl - 

(h) If —m < m < and —mto > tq, then 

*'<(!- V)4 

Proof. Let w = z^. Define h{w,fi) = 2mw — 2cw^ + 2w^hi{y/w, fi). We have w' = 
hiw^jj). Let the solution map be w{t,WQ). Since hi is a C°°, even function of z, it 
can be shown that h is C°° . Let w = ^ and we = ^e = + 0{m'^). It is easy to 

see that < if ^ is small and w < Therefore, using Taylor's formula with 
remainder, we have 



d fh{w,fi)\ _£h-w-h 



< 0. 



dw \ w J w"^ 
Here we have used the fact that h{0, fi) = 0. Similarly, since h{wE,lj) = 

d ( h{w, fi) 



dw \w — We 



< 0. 



Consider case (a), m > first, 
(i) If Wq > We, then w(t) = w(t,WQ) > we for all t > 0. Since h{wE,l^) = ^^^d 



21 



hiw, u) d , , , ^/ 9x 

^ ' < T^hiwE, /i) = -2m + O(m^) < -m, 



w — We dw 

if < m < m. Let e~^° = 1 — r^. From (w — we)' < —m{w — we), we have 

< {l-ri){wo-WE)- 

(ii) Observe that w < we- If < wq) < w, then < w{t) < w, for < t < ^. 
Since — is monotonic, 

■f/; ~ w 

= 2m - 2cti) + O(w^) 

2m ^, o. 
= 2m- — + 0(m2) 

> m, 

if m is small. Let 1 — r] — e""^"/^. Prom w' > mw, we have 

Therefore wq < {1 — ri)w(to). 

(iii) If w < w{^) < We, then for ^ < t < to, w < w{t) < we- Observe that 
^ — 3 + 0{m). Using the monotonicity olh/w, 



w 



h{w,ii) ^ h{w,n) 
w — We ~ w — We 



(2m - 2cw + Oiw^)) 

w - We) 

1 Am 



-2 + 0{my 3 
2m ^. 2x 



— + 0K)) 



m 



if m is small. Let 1 —r] — e "^"Z^. Prom [w — we)' > — yl'"^ — we), we have 

to 
2 



< 1^0 - WE\e'^"'^ 

< \wQ-WE\{l-r]). 
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Case (b), m < can be proved similar to case (a), (i). □ 

Lemma 4.4. Assume that ad — be < 0, rfi > 0,^2 > and {a + 7r'^di){d + 7r'^d2) > be. 
Then /(O = {a + ^7f^di){d + ^7T^d2)-be satisfies f (0 > for all ^ > 1 and /(O > 
for all ^ > I. 



Proof. The assumption implies that /(I) > 0. It is easy to verify that 

A 



> |{(a + ^'K'^di){d + ^TT^rfz) - be} 

From this the desired result follows. □ 

The following lemma relates the Hypothesis H4) with the breaking of the homoclinic 
orbit q{t). It is a variation of a well known result on the homoclinic bifurcation using 
Melnikov's integral. See |T6 |. 



Lemma 4.5. Consider the ODE system Let be a cross section intersecting 

the orbit of q{t) transversely. Let g(0) G S. Assume Tq(t)W^{E) r\Tq(^t)W'^{E) is one 

dimensional - spanned by q{t). Let ti > and v ± {Tq^ti)W^{E) + Tg(ti)W-''{E)} . 
Then for each k ^ koo, there exist a unique g{k) G M and a piecewise smooth solution 
U{t,k) of ([Tip that is in (-00, ti) U (ti,oo). Moreover, U{0,k) G S n WiE) 
and U{tt, k) G W'{E) with U{tt, k) - U{t^, k) = g{k)v. Here Uit^, k) and U{tt, k) 
denote the left and right limit at ti. Finally, if H^^) is valid, then -^g{k) 7^ 0. 
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5. Bifurcation equations for simple and symmetric double periodic 

solutions. 

Let X > be a small constant and S = {x = x} be a cross section that intersects 
the orbit of q{t) transversely at {x, 0, 0) G O. Assume that g(0) G S. Trajectories 
near the homoclinic orbit must hit EDO at least once. We can make S smaller so that 
trajectories starting from S must reenter O after a fixed time ti. The cross section 
S is used to fix the phase. We are not construcing a Poincare mapping: S ^ S. 

First consider a simple periodic solution of period T = to + ti. Since ti is fixed, 
the period T is determined by to- The solution can be divided into an outer solution 
U^(t) = {x^(t),y^{t),z^{t)), <t <ti and an inner solution f/* (t) = (x*(t), y*{t), z*(t)), 
< t < to- In the sequel, we use superscript (subscript) to denote inner (outer) solu- 
tions. Let the outer solution be specified by an initial value problem with the initial 
value f/*(0) = (x, yi, Zi) G S and let the solution be denoted by f/*(t; x, yi, Zi, fi). Let 
the inner solution be specified by the boundary value problem as in Lemma |4.1| with 
the boundary conditions x*(to) = x, y*{0) = yo and z*{0) = zq, and stays in O for 
all t G [0,to]- See Figure 5.1. By Lemma ^TTj such inner solution is unique and is 
denoted by 

x*{t) = x^{t;to,x,yo,Zo,lJ,), 
y*{t) = y^{t;to,x,yo,Zo,n), 
z*{t) = ^{t;zo,iJ^) + z^{t;to,x,yo,Zo,iJ^). 

Define 

x*(to,Z/o,2o,/i) = x*{0), 
y*{to,yo,zo,iJ,) = y*{to), 
z*ito,yo,zo,ij,) = z*{to), 
x{yi,zi,fi) = x^(ti;x,yi,zi,fi), 
y{yi,zi,fx) = y*(ti;x,?/i,zi,/i), 
z{yi,zi,n) = z.^{ti;x,yi,zi,jj). 



Fig. 5.1 

The end points of outer and inner solutions must match. We have the following 
equations 

Gi = x{yi,zi,ij) - x*{to,yo,Zo,n) = 0, (5.1) 



24 

yi = y*ito,yo,zo,fi), (5.2) 
zi = z*{to,yo,zo,iJ,), (5.3) 
yo = yivuzufi), (5.4) 
Zo = z{yi,zi,n). (5.5) 

See Figure 5.1. 

Substituting ( |5.4]) and (|5.5|) into ( |5.2| ), we have 

2/i = y*(^o, yiVuZuiJ,), z{yi,zuij,),fi). (5.6) 

Using the smallness of dy*/dyo and dy*/dzQ, see Lemma we solve ?/i from (|0| ) 
by a contraction principle to yield 

yi = y{to,zuiJ,). (5.7) 



Lemma 5.1. i) There exists a constant ai > such that \y\ + ^ = 0{e 

ii) y{to,-zi,iJ,) = Ry{to,zi,iJ,). 

iii) When zi = 0, (x,y{tQ,0, fi),0) E Xq. 

Proof, i) follows from Lemma 

Since {x^,(t),y^{t), z^{t)) satisfies initial values (x*(0), y*(0), z^,(0)) = (x,yi,zi), then 
(i?x*(t), Ry^:{t), Rz^{t)) satisfies initial values {Rx, Ryi, Rzi). Therefore 

w^{ti; Rx, Ryi, Rzi, fi) = Rw^{ti; x, yi, Zi,fi), 

where = x*, or z^,. Next, since {Rx*(t), Ry*(t), Rz*{t)) satisfies boundary values 
Rx*{tQ) =Rx = x, Ry*{0) = Ry°, Rz*{0) = Rz^, 

w*ito, Ry\ Rz^, /i) = Rw*ito, y\ z\ /i), 

where w* = x*, y* or z*. Based on these facts, using the uniqueness of the fixed point, 
(ii) can be verified from ( [5.6| ). 

When zi = 0, in (p. 61), let yi G Xq. Since Xq is invariant under the flow, then 
y* G Xq and G Aq, i.e., = 0. Since we can solve the boundary value problem, 
as described in Lemma [4.1| , in Xq, therefore, the right hand side of ( |5.6| ), i.e., y* is in 
Xq. We then can solve ( |5.6| ) by the contraction principle in Xq fl F . This implies that 
the unique solution y{to, 0, fi) G A^q. iii) then follows from Theorem □ 

We now substitute (|5.7| ) into ( |5.1D . Recall that /i = (£, m, fc). Gi is now a function 
of (to, ^, fc, ^i)- 

Gi(to, ^, fc, zi) = x{y{to, zi, /i), zi, /i) - x*(to, ^, /^) (5.8) 
where the arguments of y and z are (^(to, Zi, zi, ft). 
Lemma 5.2. G'i()f:o, ^, fc, — -^i) = ^1(^0, ^, ^, -^i)- 
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Proof, i) The functions y,y, z,x and x* are all invariant under the reflection R. There- 
fore Giito, i, m, k, Rzi) = RGiito, i, m, k, zi). Assertion i) then follows from the facts 
Rzi = —zi and RGi = Gi. 

ii) Set to = +00, Zi = and fi = fio where /xq = (^o? "'^O; ^00) with ttiq = (or 



(£0, ""^o) ^ r). We then have x* = from Lemma |4.1| and y = from Lemma |5.1| . We 
now show that the function Gi{oo, io, 0, k, 0) is the Melnikov function as in Lemma |175| . 

Since Aq is invariant under systems (2.1), and {x, 0, 0) e Aq, we have x*(i:i, x, 0, 0, fj,) 
e Xq. In Xq, the equilibrium E is hyperbohc with Wil^{E) = {x = 0}, Wl^^{E) = 
{y = 0}. Thus (x,0,0) G M/;"^(-E'). Consequently, U(ti), with the initila condition 
(x, 0, 0) is in W^{E) fl Aq. Observe that (1, 0, 0) is a vector orthogonal to q(ti) where 
ti is a large constant such that q(ti) has reentered O. Thus, (71(00, £0, 0, /c, 0) is 



the function g{k) in Lemma From Lemma 4^ and Hypothesis H4), we have 
^(^1(00, £0, 0, fcoo, 0) 7^ 0. Observe that Gi is a function in a neighborhood of 
(00, £0, 0, koo, 0). Thus, ^Gi 7^ for (to, ^, m, k, zi) near (cx), io, 0, fcoo, 0). □ 

Since Gi(oo, £0, 0, fcoo, 0) = 0, reflecting the existence of the homoclinic solution 
q{t) at k = koo, we can use Lemma 5^, ii) to solve k = k*{to,i,m, zi) from ( ^.8|) , if 
to ~ 00, £ ^ io, m ^ and zi ^ 0. From Lemma |5.2| , i), 

A;*(to, ^, m, —2:1) = A;* (to, £, m, zi). (5.9) 

We now substitute k = k*{to,i,m, zi) into (p. 3]), to obtain a bifurcation function 



G2(to, i, m, zi) = z*{to, y{y{to, zi,^), zi, /i), z(y(to, zi, /i), zi, yu), /i), 

(5.10) 

where /x = (£, m, /c*(to, £, m, zi)). A solution of the equation 

= G2(to,^,m,2;i) (5.11) 
corresponds to a simple period T = to + ti solution to (2.1). 

Lemma 5.3. G2{to,i,fn, —zi) = —G2{to,i,fn,zi). In particular G2{to,i,fn,0) = 0. 
The solution corresponds to zi = 0, k = k*{to,i,m,0) is a period T = to + ti SH 
solution. 

Proof. Since the functions y, y, z and z*, in the definition of G2 are all invariant 
under the symmetry R, so is G2. since G2 E Z , we have RG2 = —G2- This proves 
that G2 is an odd function of Zi. 

As in the proof of Lemma |5.2| , y(to,0,/i) G Xq. The outer solution with initial 
condition (x, y, 0) G Xq must be in Ao. Therefore the periodic solution corresponding 
to 2:1 = is in Xo- □ 

Next, we consider a symmetric double periodic solution f/(t) of period 2T. The 
bifurcation equation for the existence of such solution can be derived much the same 
way as for the simple period T solution. Therefore we only discuss it briefly. From 
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our definition, U{t + T) = RU{t),t G M. Assuming that f/(0) G S = {x = x}, we 
define 

U^t) = U{t), 0<t<ti, 
U*{t) = ult + ti), 0<t<to, 
T = to + t,. 
The matching conditions on the outer and inner solutions are 

f/*(0) = f/,(ti), 

f/*(to) = RUM- 

As before, let the outer solution f/*(t) be determined by the initial value t/*(0) = 
(x,yi,zi) and the inner solution U*{t) be determined by the boundary condition 
{x*{to),y*{0),z*{0)) = {x,yo,zo). Then we still have (U), O and Q, but (U) 
and ( ^l3|) change to 

yi = Ry*ito, yo, zq, /x), dU)' 

zi = i?z*(to,2/o, ^0,/^)- (1^)' 
Notice that the functions Ry* and Rz* have similar smallness properties like y* and 
2*. We will use the same notations as we worked on simple periodic solutions if no 



confusion should occur. As before, we can solve yi to get (|5.7|) , and Lemma |5Tl| is still 
valid. Here and afterward we use the same notations for functions y, Gi, k*, G2 when 
deriving bifurcation equations for both simple and symmetric double periodic solu- 
tions. Define Gi as in ( p.8| ). Again we have Lemma p.2| . Solving k = A;* (to, ^, 'm^, Zi) as 
before, we still have (|5.9|). Substituting k = k* into (|5.3|)^ and defining the function 
G2(to,i,'m, Zi) as in ( |5.1U| ), we have found that the solutions of the equation 

-zi = G2{toJ,m,z,) (|T3y 

correspond to symmetric double period 2T solutions to (2.1). 
Analogous to Lemma ^.3|, we have 

G2(to,i,m,-zi) = -G2(to,i,m,zi). 

However, when zi = 0, we really have obtained a SH simple period T solution tracing 
its orbit twice, since f/(t + T) = RU (t) = U{t) in this case. 
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6. Proof of the main results. 

Proof of Theorem |3.1| . Since A11A12 = (a + 7c'^di){d + TT^d2) — be, which is negative 
(positive or zero) in Q^{Q_ or F), assertions (i), (ii) and (iii) follow from the fact 
ReAi2 < 0. 

Let ((ii, ^2) e r U From Lemma WM, A„iA„2 = (a + n^7r^c/i)((i + n^7r^(i2) — be 



is an increasing function for n > 1. It follows that A„iA„2 > for n > 2. From 
ReA„2 < we have ReA^i < 0, n > 2. This proves (iv) for n > 2. The proof for the 
cases n = 0,1 are obvious and will be omitted. 

Let a = a + vr^rfi and (3 = d + 7r^(i2. a + f3 > since a + d > 0. 



Aii,Ai2 2 

^An = \{-l + {a- - (3f + 46c]-i }. 

When {di, ^2) G F, we have a(3—he = 0. Therefore ^An = \ {—1 + {a — (3) / {a + (3)} 
= -(3/{a + (3). Thus, ^ < if /3 > 0. Similarly ^ < if a > 0. It is impossible 
to have both a < and /3 < since a + /5 > 0. This proves (v). □ 

As mentioned earlier, we use only one set of notations y, k*, Gi and G2 for functions 
employed when deriving bifurcation functions for both simple and symmetric double 
periodic solutions. This allows us to treat both problems simultaneously. 

Since the proof of the main results are technical, it may be useful to preview the 
main idea used here. Consider finding a simple SN periodic solution. Recall that the 
bifurcation equation 

zi = z*{to,y,z,ij,), 

where y and z are as in ( ^.10|) , has a trivial solution = 0, which corresponds to a 



SH solution. Since we are not interested in such solution, it is reasonable to look for 
solutions of the equation z*/zi = 1. 
Let /i = (£, m, k*(to, i, ■m,zi)) and let 

zo = z{y{to, zi, /i), zi, n), (6.1) 

Then zq = if Zi = 0. We look for solutions of 

— = 1 

zq Zi 

We can show, in the limiting case, z* / Zq ~ ^/ Zq, the latter is the rate of expansion on 
the local center manifold. Also Zq/ zi k, e*{t), the latter is the rate of expansion along 
the outer solution, and its sign represents the twistedness of the homoclinic orbit. 
In particular, based on Lemma |4.2| , we can show that ^ is monotonic for zq > 
or < (Lemma |6.2| ). The proof of our main results would be easier if the outer 
solution were the multiplication by c* and the inner solution were the expansion by 
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the rate $/-2o- However, such approximations have some small errors. Care must be 
exercised to ensure that the error terms do not disturb the main terms. 
Denote 



ri{to,i,m) = I hm 2;o/zi| . (6.2) 



Then 



ri{to,i,m) = — — (to,0,/i) H h 



dyi dzi dzi dk dzi 

Assuming now to = +00, we have y = and ^ = (Lemma |5.1| ). Also 
k*{oo, io, 0, 0) = koo, reflecting the existence of the homoclinic solution q(t) at k = k^. 
From ( |5.9| ), we can show that ^A;*(oo, £, m, 0) = 0. We have shown the first part of 
the the following 



Lemma 6.1. ri(oo,£o,0) = ^(0, 0, (4, 0, fcoo)) = c*(£o), where c*(£o) is defined v 



in 



Definition p7 

Lemma |6]l| offers an easy way to compute ri(cxD, £0, 0) since c*(£o) can be obtained 



by computing the ODE system (|3.5| ) in Xi. The proof of the second equality of 
Lemma |6.1j is deferred to §7. 

In the first part of this section we assume that c*(£o) 7^ for an [q G /. Then, 
based on Lemma |6Tl|, ff^ 7^ when 2:1 = 0. We can solve Zi form ( |6.1] ) to obtain the 
inverse function 

zi = zi{to,i,m,zo). (6.3) 

Here zi is a smooth function defined for to ~ +00, m ~ 0, i ^ Eq. Assume that the 
domain of Zi is so small that 

sign{zi/2:o \ zq 0} = sign c*(£o). 

Let 

r /, . N 1- z*{to,y{y{to,Zi,fi),Zi,^),Zo,fi) , . 

Li[to,t,m) = hm , (6.4) 

20 -»o Zo 



where zi is given in ( |6.3| ). From Lemma |4.1| , we have z*{- ■ ■) = $(to? Zq, fi) + z {■ ■ 
where ■ ■ ■ represents the variables from the r.h.s. of (|6.4| ), and + = 0(6""^*"). 
It follows that 

L,(to,£,m) = ^%M + o(e--'o), a, > 0. 

OZo 

Since ^{t,Zo,fi) satisfies the equation z' = mz — cz^+ h.o.t., see (|4.3|), we have 
^$(to,0,/i) =e™*o. Therefore, 

Li(to,^,m) = e'"*o + 0(e-"i*o). (6.5) 
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Recall the monotonicity of ^{to, zq, fi) / zq, 2:0 7^ 0, proved in Lemma ^l2| . Using the 
smallness of = z*(. . . ) — ^(to, Zq, jj), we have the following results. 

Lemma 6.2. Let ±ze be the nonzero equilibria of ( H-dj ) ifm > 0. Let fiQ = {i^, rriQ, koo) 
where io & I and niQ = 0. For each < r] < 1 there exist e = e{r]) > and m > 
such that ifO < \zo\ < e, — ij,q\ < e, to > 7, and m <fn, then we have the following 
result. 

[z*{tQ,y{y{tQ,zi,ij),zi,ij),ZQ,ij)/zo] (6.6) 



dZn 



^0 



ifm > 0, 2I ^ zl, zl < (1 -r7)<l>2 



or \z. 



C6(e-™'*« - l)/zo, 

$2) $ 



$2| < {l-r])\zj 



,2\. 



ifm>0, Zq 
where — m' = lim 



z\ and e < 1 — r], 

2-^2B ,92 V /I 



, m-czl + Zohi{z^,iJ,) Zq 



2 ' 



ifm<0, and ^'^ < {1 — ri)zQ. 



Here CsjCe > \ and C7 > | are functions of (to, i,m, zq) . In all the three cases 
''9n {^(S)} = - stgn {zo}. 

Proof. Assume that rh is small so that \ze\ < e/2. 



d 

dzo 



< 



1 

W\ 



{dz^/dzo)zo\ 



d 



z\zo)-—z\0)-Zo\ + 
dzQ 

< Csup II— ^11 ■ |2;o|. 

dZn 



zo--^z\zo)- Zo\ 

CLZq CLZq 



Here the fact that |2o| < e is used. We now use Lemma ^Al We may assume that P in 
that lemma is arbitrary small at the cost of selecting a smaller e. Let < a < ao — 3/3, 
we have, 

|A(f!)|<C8|.o|e-"*°. 
azQ zq 

(i) Suppose m > 0, z^ Zq and Zq<{1 — 77)$^. Let C2 be the constant as in 
Lemma 14.21, case fii) . Then 



C2{zi - $ 



> 



C2?7$ 



>SCs\zo\e 



-ato 



if to > 7 and e > are sufficiently small. Here we have used the facts $^ > Zq and 
small. Therefore \:f-{-)\ < \\^{-)\- Since C2 ^ 1, C5 > i. Therefore, 



IZjp Zn 



IS 



p^a) follows from Lemma case (ii) 
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Suppose now m > 0, z% ^ Zq and \z% — < (1 — ri)\z^ — Zq\. Then 

> 1 - (1 -r^) =7]. (6.7) 

|$| 
C2V— 

'4 

> 3C8|zo|e-°*° 

The last inequahty is based on the fact |$| > l^ol and \zo is small. Based on a similar 
argument, ( |6.6| a) follows from Lemma [4.2| , case (ii). 




(ii) If m > 0, Zq = z\ and e < 1 — r], then 

|(e-'"'*°-l)/2o| > Wl^ol >3C8|^o|e-"*'' 



if |2;o| < e and to > e > is sufficiently small. As before, ( |6.6| b) then follows from 
Lemma |4.^, (iii). 

rj)zQ, then 

Ci(22 _ $2) $ 



(iii) If m < and $2 < (1 



> 



m-czl + z^hi{zo,fi) zl 



m- cz^ + ZQhi{zo, n) 
> 3C8|zo|e-"*«, 

When deriving the last inequality, we assume e is sufficiently small so that sup^{h{z)} 
= rh with m < a. Then |$| > |2;o|e~'"*'' by the Gronwall inequality. Therefore, if to is 
sufficiently large and m — cZq + ZqUi^zq, 11) is small, the last inequahty holds. (|6.6| c) 



then follows from Lemma 4^, case (i). 



The proof of Lemma 6.2 has been completed 



□ 



Lemma 6.3. Under the same conditions of Lemma \6.2^ , in any of the three cases, 
dfj), dp ) or (\6. Of ), if we choose smaller e > 0, we have, (see \6.d[ ) for Zi), 

d fz*\ I >0, 
<0, 




if ~zi < 0, 
if Zl > 0. 



d fz*\ d 
Proof. Observe that 7- — I — ) = -7: — 

OZq \ZiJ OZq 

is an odd function of Zq, we have 



We first show that 



d 




dzQ 


(1)1 


d 




dzo 





£2£o 
^z^ 



Zl dzo \Zq/ Zq dzo \ Zi 



z* d fzo', 

— ■ omce Zl 



[Zl-T^Zo)/zl\ <Cg\z^o\/^f- 
OZq 



> 2Co 



Zl 



since this will imply that the sign of 



a 

dzo 



(^j^^ is determined by the sign of (f^)' Based on |$| > \zo\e and \z^\ < 
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C\zo\e 



-aito 



we have for large to, \z*\ < 2|$|. We then need to show 



_d_ 



(6.8) 



where Cio = 4C9C11 with Cu = sup 
(i) li m > 0, zl; ^ Zq and Zq<{1 — rj)^'^, then from (I 



I zi I 



_d_ 

dzo \zo 



|r/$| 



2 kl; z'qI 



provided that \z% — Zq \ < which can be achieved by choosing smaller e. 
(ii) If m > 0, z| 7^ z'^ and {z^ - < {1 - r])\zl; - ^^\, then from (|6l6H) and 

\ri<!?\ 



d 



> 



2zl 



> Cio\<l> 



dzo \zo, 

provided that z^ < 20^1 which is valid if e is sufficiently small, 
(iii) If m > 0, 2;q = z% and e"*"'*" < 1 — 77, then from (|0|b), 



d_ 

dzn 



V^o 



> 



2\zo\ 



Here we need \zq^\ < 



2C 



which is valid if e is small. 
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(iv) If m < and $^ < (1 — 77)2:0, then from ( |6.6| 



d 



z" 



dzo \zo 
2 



> 



\7]<l>\ 



> Cio|$ 



2\m - czl + 4^1(^0, At) I 
provided that \m — czq + . . . |< that is valid if e is small. 

In all the cases, the sign of ^ (^C^ agrees with that of (f^) • From Lemma |6.2| 



sign 



d 



dzo \zo 

Therefore sign{^ = - sign{zi} 



sign {zo}. 



□ 



Corollary 6.4. Under the conditions of Lemma 6.i, we have I — is a decreasing 
function of \zi\. 



Recall that ti is fixed and T = + depends solely on t^. Let 

L{T,i,m) = Li{to,e,m)e"''\ 
r{TJ,m) =ri(to,^,m)e-""*^ 



(6.9) 



We shall prove that with such functions L and r. Theorem 3.4 is valid. 
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If no arguments are given, for notational simplicity, $ means $(to, -Zo? /^); is 
defined in ( |6.1| ), z* = G2(tQ, £,171, zi) is defined in ( |5.10| ), zi = z\ is defined in ( |6.3| ), 
^ means z^{to,y{y{to,zi, fi),zi, fi), z'^, fi) = 2;* — $. Assumptions following a case 
number are valid until a new case is encountered. 



Proof of Theorem 3.4. From 



if m = 0, 



r(cx),£o,0) = ri(oo,£o,0). 

Thus, Theorem |3.4| , 1) follows from Lemma |6.1| . 

From ( |6.9| ) again, it is easy to see that 2) follows from (|6.5[). 
We now prove Theorem |3.4] , 3). From ( |6.4| ) and Lemma |4.1| , 



Li(to,i,m) 



+ 



we have 



dzodm 
«9^fa,0,/i) 
(9^0 

dLi{to,i, m) 



dzo 



9zn 



ai > 2/? > 0, 



dm 



toe"''° + 0{e 



(-al+2/3)^o^ 



Now that 



d 



„ -(Liri) = f^ri + and < C, we have 7r-(^in) = toe'^ 

am / ym ^ ^ am ' am ' — ' am \ ^ / " 

ri(to,^,m) + 0(6"^*°). Let O = {|£ - 4| < 6, \m\ < 6} and let to > 1/6. If 5 > is 
sufficiently small, we have ri{tQ,i,m) = Cic*(£o) where Ci > ^. Thus, 



sign 



_d_ 

dm 



sign {c*(£o)} 



The assertion in 3) follows by observing that L ■ r = Li ■ ri. 

Consider c*(£o) > and look for SN simple periodic solutions first. We only need 
to solve ( |5.11|) for zi > Q since G2(to,i,m, Zi) is odd in zi. If zi is sufficiently small 
and to is sufficiently large, we have ri{to,i,m) > and zq/zi > 0. Thus we need to 
consider 2:0 > only. We now look for a solution zi G (0, Q with the corresponding 
Zq G (0, e) where (— e, e) is the coordinate chart in the 2;-axis for Wfg^{E). 

Since Zq = if zi = and zq depends continuously on zi, we can choose a small 
constant C > so that zi = ( implies zq < e. If 6, which defines the set O, is small, 
then either 2;^; < ^/2, m > or does not exist (m < 0). We can choose t > 
so large that > t implies that < $ < 3^/4, and \z'^\ < C/4. The first estimate 
is based on ^ ze < C/2 or as to oo- The second estimate uses Lemma ^A . 
Therefore (^2(^0, m, C) = -2* = $ + -2* < C- 
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Suppose now Li{to,i,m)ri{to, £,171) > 1. Then 



z* 



lim — — - > 1. 

20-0+ z^ z^ 

From this, there exists a small < Zi < ( such that G2(to, i, m,zi) > z\. Thus, there 
exists at least one solution < 2:1 < C for (|5.11| ) if Lr = LiVi > 1. 

In the rest of the proof, we discuss the uniqueness or nonexistence of SN simple 
periodic solutions. Please refer to Figure 4.1 for the flow on Wf^^{E). 
Case (i): < c*{io) < 1. For any < r^i < 1, by choosing smaller 6 and e, 

0<Vi<4/4<^-Vi- (6.10) 

Let m < and zq > 0. Then < ^ < zq and \z^\ < 2oe""^*°. If to is sufficiently 
large, from Lemma [4.1|, we have {z*y/zQ < 1/(1 — 771). Combine this with (|6.10|) , we 



have {z*y < zf. Therefore, there is no solution for (|5.11|) . 

Let m > and zq > ze- Similar to the previous case, we find no solution for (|5.11|) . 
For m > 0, define 



'1 

.2 



Zm = SUp{zo \ Zq< {1- T])^ ,0< Zo< Ze} 

for some < rj < rji. Clearly < zm < ze- 

Let m > and zq G {zm^ ze)- We have Zq>{1 — 7/)$^. By choosing a smaller 5 
we have Zq > {1 ~ r]i){z*y. From ( |6.10|) , we have (-2*)^ < zf. There is no solution for 



(|5.11| ) in this case. 



Let m > 0, and zq G (0, zm] if zm > 0. At zq = zm we have ^ (^^^ < 0, based on 
Lemma [4.2| . We infer that Zq < {1 — 77)$^ for all zq G (0, zm]- Lemma |6.3| implies that 
{z*/zi) is strictly decreasing in that interval. This proves the uniqueness of solutions 
of ( ^.llj ) in this case. 

Let m > and 1 > Li(tQ, i,m)ri{tQ, i,m) = \imz-^^^Q+{z* / Zi). In the case Zq G 
(0, Zm], arguing as in the previous case, we have z* /zi < 1. This proves the nonexis- 
tence of solutions of (|5.11| ) if Lr < 1. 



Case (ii): C*{io) > 1. For any < r^i < 1, by choosing smaller e and 5, 

zf<il-ViK (6.11) 

Let m > and < Zq < ze- Since ^ > zq and \z^\ < ZQe~"'^*°, we have {z*Y > 
(1 — ?7i)zg if to is large enough. From (|6.11|) , there is no solution to (|5.11| ) in this case. 



For m > 0, define 

Zm = inf{2;o | "^"^ < (1 — f])^o ^^^^ ze < zq < e] 

for some < r/ < r/i. From the phase diagram, cf. Figure 4.1, 2;^; < 2;^ < e if r/ is 
sufficiently small. 

Let m > and Zq G {ze^z^)- Then $^ > {l-7])zl. We can have {z*Y > (1-^71)^0 
if we choose to large enough. Then z* > zi based on ( |6.11|) . Equation ( ^.11|) has no 
solution in this case. 



34 



Let m > and zq = Zm- Then $^ < (1 — rfjZQ. This imphes that — zW < 
{1 — Tj) \ zl — zW- We then can show that / Zq is decreasing for zq G [zm,e]- Thus 
$^ < (1 — rj)zQ for 2;o G [^m^e]- From Lemma z* / Zi is strictly decreasing. The 
solution to ( |5.11| ) is either unique or does not exist. We show Li(to) TT^)^i{to, i, m) > 
1 in this case, so that the nonexistence becomes impossible. In fact, if e > is small, 
then riitQ.l.m) > 1 + 772 for some 772 > 0, due to c*(£o) > 1- Because m > and 
Li{to,i,m) = e™**' + O(e""i*o), let to be sufficiently large, we have Li > (1 + 772)"^ 
Therefore Liri > 1. 

For m < 0, define 



'5 

z"" = inf {zo I < (1 - < zo < e} 

for some rj < rji and 77 < 1/4. Clearly < < e if 77 is sufficiently small. 

Let m < and 2;o ^ (0,2"^) if 2;"^ > 0. Then <l>^ > (1 - r])zQ. We can make 



\2 



> (1 — r7i)2Q by choosing e < smaller. Therefore Zi < z* by (|6.11 ). There is 



no solution to (|5.11|) . 

Let m < and Zq = z™. Then $^ < (1 — ri)zQ. By Lemma [4.2| , ^^/-Zq is decreasing 
for Zq G [2;"^, e). Thus $ < (1 — ri)zQ for 2:0 G [-z™, e]. Then is strictly decreasing 

by Lemma |6.3| . There is either no solution or the solution is unique to ( |5.11| ) when 

Zoe[z"',e). 

Let m < and < Li{toJ,m)n{toJ,m) < 1. By (l6Tl|) , rj > {1 - 7]i)-\ Thus, 
Li < 1 —rji. If we choose e > small, we have {z*y/zQ < 1 — 772 for some < 772 < Vi- 
And also ^"^/zq < 1 — 77 for some < 77 < 772 in the interval [z'",e]. Thus z*/zi is 
strictly decreasing. Since lim^-^^o ^ there is no solution to ( 5.11 ) in this case. 
Case (iii): C*(£o) = 1- For any < 77 < 1, we can choose a smaller e so that 

l-r]<{zo/zi) <l + r]. (6.12) 

If e is small, then ri(tQ,i,m) < 1 + 77. Let Li(tQ, i,m)ri(tQ, i,m) > 1 + 5 for some 
6 > f]. Then 

1 + '^ ^ . 

Ll > —— > 1 + 7/1 
1+7/ 

for some < rji < 6. From Li(to,i,m) = e™*° + 0(e"°^*°), if to is large enough, we 



have 777, > 0, and 777,^0 > ei for some ei > 0. From Lemma [4.3| , case (a), we have 
Zq < {1 — 773)$^, or 1$^ — < (1 — t]3)\zq — z%\ or e™*° — 1 > 773 for some 773 > 0. 



Therefore z* / zi is strictly decreasing and the solution to ( |5.11|) is unique 



Let Li{to, i^m)ri{to, i,m) < 1 — 6, for some 6 > rj. By a similar argument, 7?7-to < 
— ci for some ei > 0. Also we must have 7ti < 0. From Lemma [4.3| , case (b), we then 
have $^ < (1 — 7/3)^0 for some 7/3 > 0. Thus z*/zi is strictly decreasing. There is no 
solution to ( |5.11|) since lim2j^o+ z* / zi = Li ■ ri < 1 — 6 . 

We have completed the discussion for the case C*{io) > 0. 

Consider C*{io) < and symmetric double periodic SN solutions next. We need to 
solve (^[11])'. We can divide the case into three subcases — case (iv): —1 < C*(£o) < 0, 



case (v): C*{io) < —1 and case (vi): C*{io) = —1. They are analogous to cases (i) 
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(ii) and (iii) respectively. The proofs are similar to the previous cases and will not be 
rendered here. 

This completes the proof of Theorem p.4| . □ 

Proof of Theorem |3.5| . From our assumption, c*(£o) = ''"1(00, ''^o) = and 
-^ri(oo,io,mo) 7^ 0. Using the imphcit function theorem we can find a unique 
function i = 



t{m) so that r(oo, £*(m), m) 
Since |y(to, ^i,/^) |< Ce~"i*" and |fc*(to, 



0, \m\ < S. 

m, Zl) — k*{oo, 



|< Ce 



we 



0{e 



Please refer to (|6.1|), (|6.2|) for the definitions of Zq 



have ri{to,£*{m),m 
and ri. 

Since zq is an odd function of zi, for some C* > 0, we have 

^° <C{e-''''' + zt+\i-t{m)\). 



Since ^(to,zo,fi) satisfies the equation z' 



\zo\ < e and |$| < e. Thus, from Lemma 

z 



mz — cz^ + ■ ■ ■ , we have 1 — 1 < e™'*" if 



Zo 



if to > ^ is sufficiently large. We now choose S{T) = Ce "^'^ where C is a small 
constant. U\i-t{m)\ < 5{T) and Zi < {5{T)Y/^, then C{e-"^'° + zl + \i- t{m)\) < 

< 1. The bifurcation equation (|5.11| ) or (|5.11|) ^ has no 



ig mto^ Therefore, 
solution in this case 



si 










Zl 



□ 
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7. Numerical test on a predator-prey model 

The following predator-prey model was proposed by Freedman and Wolkowicz [13] 
to describe group defense of prey against predatation. 

ii = 2u(l — — ) — 9vp(u), 
k 

V = v{--f + 11.3p{u)), (7-1) 
u = prey, v = predator. 

where p{u) = u/{u'^ + S.SSm + 13.5) represents the interaction between prey and 
predator. For a large range of ipf^k), (|7. 1| ) has two interior equilibria {uq^Vq) and 
{uq,Vq), with Mo < Uq. Here p{uq) = p{uo) = while Vq, Vq can be solved from 
the first equation of (|7. 1|) . We are interested in the equilibrium {uo,vo), which is 
hyperbolic, and shall be denoted hj E = E{'j, k). 

Let the Jacobian matrix at be ^ ^ ^ ^ . It is shown in [25] that a > 0, 6 > 

0, c > and d = 0. Thus, Hypothesis H^) in §3 is satisfied. 

Freedman and Wolkowicz have discovered a curve 5 C such that if (7, k) G S, 
then (|7. 1| ) possesses a homoclinic solution q{t) asymptotic to the equilibrium E{'y, k). 
Numerical computation shows that the curve S can be parameterized by uq and is 
plotted in Figure 7.1. For each (7, A;) G iS, let 7 be fixed and let k vary. Then the 
homoclinic solution breaks. The derivative of the gap between W^{E) and ^^{E) 
with respect to k can be evaluated by the Melnikov integral M^{k), as in H4). The 
Melnikov integral has been computed numerically, and the result is plotted in Figure 
7.2. Evidently, M > for all the values considered. Thus, Hypotheses H2) and H4) 
in §3 are satisfied for those parameter values. 



Fig. 7.1 



Fig. 7.2 
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The smooth dependence of M^{k) on mq indicates that M > is not a numerical 
artifice. 

In the remaining of this section, we fix {uq, 7, k) = (5.49178,1.0,6.87433). After 
adding diffusions {diu^^, d2V^^), we consider a system of PDEs in the domain < ^ < 



1 with Neumann boundary conditions, cf. ( |1.2| ). Let T be the curve in (c/i, (i2)-plane 
on which (1.2) has a zero eigenvalue with associated eigenvectors in Xi. Since bc> 
and a > 0, r is depicted in Figure 3.1, Case 1. 

For {di, ^2) e r, we now compute Wf^^^E) and the flow on it, up to 0{p^) where p = 
\u — uo\ + \v — Vol- Since the boundary conditions are of the Neumann type, we will ex- 
pand (m, v) into Fourier cosine series. Let (mq + I]o° ""n cos nvr^, vq + f „ cos nvr^ ) G 
Wig^{E). Let (mi cosvr^, fi cosvr^), with {ui,Vi) = (77, 1), be the unique eigenvec- 
tor corresponding to the eigenvalue 0, up to a constant multiple. Then ui = rjVi + 
0(ui), Un = W^ivi) and w„ = f*(fi) for n 7^ 1, where 0, m*, f* = O(f^). This 
is due to the fact that Wf^^{E) is tangent to the zero eigenvector corresponding to 
{ui,Vi) = {rj, 1). Because of the R symmetry (Theorem 2.7), we have = 0{p^). 

The Fourier coefficients (u„, f„) are functions of t. They satisfy 



-din^T^'^Un - aUn - bVn + [/(m, V)]n + 0{p'^), ,^ ^ 

V'n = -d2n^T['^Vn - CU„ + [g{u, V)]n + 0{p'^), 



where / and g are polynomials of degree 3. For h G ^^^(0, 1), we use [h]n to denote 
the nth Fourier cosine coefficient for h. Using some basic trigonometry formulas, we 
can rewrite [/(-u, f)]„ and [g{u,v)]n in terms of {un}'^, {vn}'o'- Only finitely many 
terms are needed here since other terms will be included in Olp"^). 

We can now use the Taylor expansion method in [2] to obtain a power series ex- 
pansion of (f){vi) and the fiow on the center manifold. The function has the form 
(j){vi) = cvf + 0{vf). And the flow on the center manifold has the from 



v[ = cvl + 0{p^). 



Fig. 7.3 

When {di, ^2) moves along the curve F, values of c have been computed numerically 
and the results are depicted in Figure 7.3, with c against vr^rfi. It verifles that c < 
for the portion of F under consideration. Since there is a diffeomorphism between vi 
and z, Hypothesis i^s), in §3 has been verifled numerically. 
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Fig. 7.4 

We now compute c*{di,d2) for {di,d2) G F, < vr^rfi < 3. Again, we fix 
Mo=5. 49178, 7 = 1 and /c=6. 87433. Numerical results of c* are depicted in Fig- 
ure 7.4. We have found a point (i^=0.183 such that c*{di,d2) < (> or =) if 
di < (> or =) dl. The results also show that dc*{io)/de ^ where (4,0) G F 
corresponds to di = dl- Therefore all the twisted, nontwisted and degenerate cases 
have be found in Freedman and Wolkowicz's example. However the case c*{io) > 1 
or < — 1 has not been found in this example. Numerical and theoretical result also 
indicate that there is a point {di,d2) = (0.0093,0.0093) where c* = 0. However, 
the numerical error near that point is too large to be trustworthy. Thus, we do not 
include it in Figure 7.4. 



We end this section by proving Lemma 6.1 



Proof of Lemma 6.1. Recall the definition of z in §5 and ri in ( |6.1| ) and ( |6.2| ). 

We need to consider the z-th. component of ^f/*(ti, x, yi, zi, /i), with yi = and 
zi = 0. Let -^U^{t) = {x{t),y{t), z{t)), < t < ti. It satisfies the linear variational 
equation (3.4) and the initial conditions are x(0) = 0,y{0) = 0,^(0) = 1. We now 
extend the solution {x{t),y{t), z{t)) to t < 0. Notice that we are treating an infinite 
dimensional system, so the backwards extension of a solution is not unique. Using 
the fiat coordinates (4.1), in a neighborhood of 0, we write (3.4) as 

x' = Aix + {D^gi)x + {Dygi)y + {D^gi)z, 

y' = A2y + Dyg2{xq{t),yq{t),Zq{t),ij)y, (7.3) 

z' = A3Z + {D^g'i)x + {Dyg'i)y + {D^gs)z, 

where q(t) = {xg(t),yg(t), Zg{t)) in the fiat coordinates. Here we have used the facts 
yq(t) = 0, t < 0, and g2{x, 0, z, fi) = to simplify the second equation of ([7.31). 

First let y{t) = for t < 0, which solves the second equation. Then {x{t), z{t)) can 



be solved uniquely from (|7.3|) backward in time. We now show x{t) = for t < 0. 
If (x(t),0,z(t)) is a solution for (^), so is R{x{t),0, z{t)) = {x{t),0, -z{t)). Thus, 
(a;(t),0,0) is a solution of ( |7.3| ). Since a;(0) = 0, solving the one-dimensional ODE 
for x{t) we have x{t) = for t < 0. Observe that gi{0,y, z, fi) = 0. Thus, the first 
equation is valid even if z 7^ 0. The equation for z{t) becomes 

z' = A3Z + D^gs{xq{t),yg{t), Zq{t), ij,)z, (7.4) 

with A3 = 0. In our flat coordinates, Zq(t) = and yq(t) = 0, we have 
Dzg3{xq{t),0,0, fi) = Dz(73(0, 0, 0, /i) = 0, since zero is an eigenvalue for (^1,^2) G F. 
Here we used the fact gsi^x, 0, z, jj) = (73(0, 0, z, fi) on Wf^^. Thus, (|7.4|) becomes z' = 0, 
and z{t) = 1 for t < 0. (x(t), y{t), z{t)) = (0, 0, 1) for t < 0. 

We now have (x(0), y(0), z{0)) G Xi and shall remain in Xi for t > 0. In particular, 
x{t) = for all t G R. According to §3, {x{t),y{t), z{t)) (0,0, c*) as t +00. 
However, because the coordinates are fiat, (73(0, y, z, jj) = (73(0, 0, z, /i). Also Xq{t) = 
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for t > ti. Thus ^g3{0,yg{t),0,fi) = for t > h. Again z{t), t > h satisfies (|7j) 
with Xg{t) = and Zg{t) = 0. Since ^^^3(0, 0, /i) = D,gs{0,0,0, /j,) = 0, We 
have z{t) = constant for t > ti. Thus z(ti) = c*. This proves Lemma |6.1| . □ 
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FIGURE CAPTIONS 

Figure 1.1. Twistedness of the homoclinic orbit is determined by comparing 0(— T) 
and 0(T). 

Figure 3.1. The curve F divides the first quadrant of the {di, ^2) plane into two 
parts. All the possibilities are listed, except for permutations of di and ^2- 

Figure 3.2. A sketch of the bifurcation diagram in the {i, m) plane. SN simple or 
symmetric double periodic solutions occur in the shaded 

Figure 3.3. A sketch of all kinds of homoclinic, heteroclinic and periodic solutions. 
Figure 4.1. Phase diagrams for the flow on the center manifold. 
Figure 5.1. A sketch of the inner solution U* and outer solution U^. 

Figure 7.1. Values of (7, k) where a homoclinic orbit to (7.1) exits are plotted, 

using Mo as an independent variable. 

Figure 7.2. The Melnikov integral M^{k) is plotted, using uq as an independent 
variable. 

Figure 7.3. Values of c are plotted when [di ^2) moves along F, which are parame- 
terized by di. 

Figure 7.4. Values of c* is plotted when [di ^2) moves along F, which is parameter- 
ized by di. 



